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Abstract 

We study a stochastic multiplicative system composed of finite asynchronous ele- 
ments to describe the wealth evolution in financial markets. We find that the wealth 
fluctuations or returns of this system can be described by a walk with correlated 
step sizes obeying truncated Levy-like distribution, and the cross-correlation be- 
tween relative updated wealths is the origin of the nontrivial properties of returns, 
including the power law distribution with exponent outside the stable Levy regime 
and the long-range persistence of volatility correlations. 
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1 Introduction 



Multiplicative processes have been well studied in different contexts and widely 
applied to various research fields, such as the biological, social, and economic 
systems [1-6]. One of the major interests in these processes is the generation 
of power laws [1,3-5,7], which have been observed in many natural domains 
and indicate the properties of scale invariance and universality. In fact, the 
investigations of power law behaviors in systems with stochastic multiplicative 
dynamics can be traced back to decades ago, and the underlying physical 
mechanisms are still to be understood. 
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For the applications in financial problems, the assumption of multiplicative 
property for the individual capital investments Wi (the index i = 1, ...,N may 
correspond to various investors/ traders or companies (stocks) in the market) 
leads to a term 

Wi{t + 1) ^ X{t)wi{t) (1) 

appearing in the dynamics of wealth/capital evolution. That is, the individual 
capital Wi at time t+1 is proportional to the invested capital itself. The random 
factor \{t) reflects the relative gain/loss incurred by individuals between time 
t and t + 1, and is chosen from a probability distribution 11 (A). 

The Pareto power law distribution [8] of individual wealths w 

P{w) ~ w-^-""', (2) 

has been found in many of the previous studies of multiplicative processes 
[1,3-7], which all follow the above dynamics (1) (plus some additional crucial 
dynamical elements that we will discuss later), and been verified repeatedly in 
the last hundred years in most of the capitalistic societies, e.g., for individual 
income and wealth [8,9], size of business firms [10], etc. In spite of the signifi- 
cant fluctuations in the total wealth (with occasional spectacular booms and 
crashes), the exponent was observed typically in the range 1 < q;^^ < 2, 
i.e., within the stable Levy regime [11]. 

Besides the above distribution property of individual wealth, we are more 
interested in studying the relative fluctuations of the total wealth W{t) — 
J2iLi'Wi{t), which may represent the capitalization (total market value) of a 
company, or the market price of a stock (when normalized to the total number 
of shares that stock has on the market), or the market index. The fluctuations 
of W{t) over arbitrary periods of time r deflne the market returns r(r) during 
these time periods: 

r{T) = In W{t + T) -In W{t), (3) 

which have attracted much attention in recent flnancial studies, in particu- 
lar for their statistical properties [12,13]. There are some empirical features 
generic for different flnancial markets, including: 

• The power-law decay of probability density of returns r 

P(r) ~ r-^-", (4) 

in tail region, but with exponent a > 2 different from that of wealth distri- 
bution (2), i.e., outside the Levy regime [14-16]. For larger r values, P(r) 
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can be fitted by an exponential in some observations [17-19]; 
• The very short range (about a few minutes) correlation of returns and long- 
range persistence for correlation of square or absolute value of returns (so- 
called volatility clustering) [12,13,15,17]. 

It is interesting to not only reproduce these stylized facts through microscopic 
models, but also understand the intrinsic mechanisms dominating the pro- 
cesses. 



2 Generalized Lotka-Volterra system 

The framework we use to describe the dynamics of financial systems is the 
Generalized Lotka-Volterra (GLV) model proposed by one of us a few years 
ago [3,4]: 

Wi{t + 1) = X(t)wi(t) + aw(t) + c{wi, W2, wn, t)wi(t), (5) 

with the average wealth iju{t) — W{t)/N. The multiplicative property at the 
individual level, i.e., Eq. (1), is expressed in the first r.h.s. term, and the 
second term represents the property of wealth at the social level, which may 
correspond to the social security, subsidies, or funded services, and can also be 
interpreted as arising from the diffusion of wealth between agents (by services, 
taxes, etc.). This term is set to be proportional to the average wealth id with 
an important coupling parameter a, supplying the correlation and coupling 
between investors or companies that is crucial for the distribution properties 
of wealth and wealth fluctuation shown below. Generally, the last term cwi 
which corresponds to the competition in the market, does not qualitatively 
affect the properties of the model [20] . For simphcity, we set it in the present 
study. Therefore, we have 

Wi{t + 1) = X{t)wi{t) + aw{t). (6) 

Note that in this stochastic multiplicative system an asynchronous updating 
mechanism [3,4,21-23] is used, due to the fact that in a real market the capitals 
of different agents or companies are not updated simultaneously. Thus, at each 
discrete time one of the elements i is chosen randomly and updated according 
to Eq. (6), while the other Wj's keep unchanged. Moreover, the random factor 
A is taken in a rather narrow range around 1, since in practice the price returns 
for very small time intervals are usually rather small. 
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3 Power law behaviors and cut-off effect 

The system described by Eq. (6) does not approach a steady state; instead 
the total weahh grows (or decays) exponentially with superimposed fluctu- 
ations, i.e., W{t) ~ exp{Kt/N) with k depending on a as well as the mean 
and standard deviation of A distribution. Nevertheless, the instantaneous val- 
ues of wealth (or the relative value Wi{t)/W{t)) fulfill the Pareto power-law 
distribution (2), with exponent < 2 [21], that is, inside the stable Levy 
regime. 

More interesting results for this multiplicative system are obtained by studying 
the wealth fluctuations (i.e., returns) r(T) deflned by Eq. (3). (Note that due 
to the growth effect of total wealth W (~ exp{Kt/N)) in this model, the 
return in (3) should be detrended, that is, change by a constant: r(T) — > 
r(r) — kt/N .) For a certain time interval r, the return r(r) can be written as 
E*tr'ln[(l^(s+l)-l^(s))/W^(s) + l]=Esln[(A(s)-l)tf;i(s)/l^(s) + a/iV+l] 
according to Eq. (6) [24]. Since the relative value of wealth Wi/W at each step 
is small and generally A — 1 is in a narrow range around 0, one obtains that the 
properties of the (detrended) return r(T, i) are the resultant of r steps walk 

t+T-l 

R{r) = E iKs) - l)w,{s)/W{s), (7) 

s=t 

with step sizes governed by power law distribution (2) with cut-off effect. 

For r = 1 the distribution of returns differs from that of the relative wealth 
Wi{t)/W{t) only by a random factor \{t) — 1, and then displays a similar 
power law behavior (4) with exponent a ~ q;^„ and a sharp peak around 0. 
As T increases, the consequence of r step walks described by Eq. (7) is the 
smoothening of this tip. 

Note that for any flnite N the distribution of Wiit) /W{t) is not a perfect power 
law even if one assumes that of Wi to be so, since this distribution is truncated 
from above: Wi/W < 1, shown as a bending in a log- log plot for very large 
values of Wi [25]. We have proven [22] that this upper cut-off effect leads to 
an exponential decay for extremely large values of returns and for very wide 
range of interval r. 

Thus, the process described by Eq. (7) looks like a truncated Levy walk leading 
to similar power laws for different r and slow convergence to Gaussian [26]. 
However, our studies below indicate a more complicated situation, that is, this 
process is different from the ordinary Levy walk due to the coupling between 
step sizes. 
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Fig. 1. The probability distribution of returns r(T) defined by Eq. (3), for N = 500 
and different time intervals r from A'^ to SOOOA?^. A uniformly distributes between 0.95 
and 1.05, and a = 0.0002. The results are averaged over 1500 runs for r < 10007V 
and 600 runs for r = 5000A'".(a) Semi-log plot, where for the largest r (= 5000iV) 
Eq. (7) is used; (b) Log-log plot for positive returns. 

The numerical results of Eq. (6) for the distribution of returns r(T) (3) are 
shown in Figs. 1 and 2, where the random gain/loss factor A has an uniform 
distribution in the range 0.95 < A < 1.05 and the coupling parameter a = 
0.0002, with the corresponding power-law exponent of wealth distribution a^, 
about 1.5 as in usual financial application. (In all the numerical simulations 
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Fig. 2. Log-log plot of the probability distribution of positive returns r(r) (averaged 
over 100 runs), for system size A'^ = 5000 and different time intervals t = N, 5N, 
5m, and lOOOiV. 

of this paper the results are calculated after t = lO^N updatings, so that the 
returns series are stationary after detrending the wealth growth factor.) As 
expected according to the above analysis, the semi-log plots of distribution in 
Fig. 1 (a) exhibit the smearing out of the sharp peak for small time intervals 
into a dome-like shape for large one, and one obtains the exponential-type 
behavior for extremely large r due to the cut-off effect of finite N system, 
even for very large interval r, which has been observed in some empirical 
studies [17-19]. 

The power-law scaling region described by Eq. (4) can be found in the log-log 
plots of the distribution (see Figs. 1 (b) and 2), and interestingly, exponents 
much larger than that of the wealth distribution (2) are obtained. For small 
interval r, the exponent a is within the stable Levy regime, i.e., a^, < a < 2, 
while for large intervals (r > N) we obtain the result consistent with the 
recent empirical observations for both the stock index (German DAX [14] and 
S&P 500 [15]) and individual stocks [16], that is, the effective exponent a > 2 
outside the Levy regime. Similar to the findings of real markets [14-16], the 
exponent increases very slightly with the increase of the time interval, due to 
the very slow convergence to Gaussian, as shown in Figs. 1 and 2. 

The extension of this power-law region with exponent a > 2 depends on the 
size of the multiplicative system. By comparing the results in Fig. 1 for 
A^ = 500 with that in Fig. 2 for A^ = 5000, one can find that the range of the 
power-law scaling is longer for larger systems. This phenomenon is attributed 
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Fig. 3. Log-log plot of the probability distribution of positive returns r(T), for 
N = 500 and different a values 0.0003 and 0.0002. Other parameters are the same 
as Fig. 1. 

to the exponential truncation effect discussed above for finite A^, which appears 
more obvious for smaller system. As observed in the log-log plots, the deviation 
and bend-down from a straight line occurs for large r, which has been found 
very recently in the Hang Seng Index (HSI) of Hong Kong [17] if one skips the 
data of first 20 minutes in daily opening. The similar effect of finite size N 
has been found in some microscopic models, in particular the Cont-Bouchaud 
percolation-type model [27], where the a > 2 power-law exponent is obtained 
over an intermediate region which becomes infinitely long for infinite large 
market size and computer time, while the asymptotic return distribution is 
expected to be a (stretched) exponential decay in the tails. 

However, different from the percolation model, the exponent a for large in- 
terval r decreases with increasing size N in our system. As shown in Figs. 
1 (a) and 2, a ~ 2.7 for A^ = 500, and a ~ 2.3 for A^ = 5000. Moreover, 
for A^ = 10^ (not shown here) a ~ 2, which is the border of the stable Levy 
regime. Therefore, the property of returns distribution is size-dependent, and 
one can obtain that the effective N for real market is not very large according 
to the high empirical a value. 

Although the value of the exponent for wealth distribution depends on the 

parameters of the system, the behaviors of the a > 2 power-law scaling for 
return distribution can exist for different parameters (e.g., coupling parameter 
a and random factor A) corresponding to 1 < < 2. When increasing the 
value of a, one still obtains the a > 2 power-law behaviors with larger a value 
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and longer extension for the same A^. In Fig. 1 N = 500 and a = 0.0002 (with 
the corresponding a^y about 1.5), we have a ~ 2.7 but with a rather hmit 
range of power law regime, while for larger a = 0.0003 as well as the same N 
and range of A (with about 1.7) a — 3.1 is obtained, and more importantly, 
the power law regime is much longer (close to 3 orders of magnitude) , as shown 
in Fig. 3. 



4 Properties of correlations: Correlated walk steps and volatility 



These results of o; > 2 are nontrivial, and different from the expectation that 
the wealth Wi and the fluctuation r should have tails obeying the similar 
power law, since a random walk similar to Eq. (7) with steps of sizes domi- 
nated by Levy-like distribution with finite variance still leads to a behavior 
within the stable Levy regime before the crossover to a Gaussian process [26]. 
However, this expectation is only valid for the random walk with statistically 
independent step sizes, while for the multiplicative process (6), there are cross- 
correlations between the relative wealths Wi{s)/W{s) comprising the returns 
in Eq. (7) at each micro-step s. This can be verified by studying the properties 
of time series Wi{t) /W{t), where Wi refers to the newly updated wealth at each 
time t. As expected, the probability distribution of the walk process of Eq. (7) 
is almost indistinguishable from that of return r(r) shown in Figs. 1,2, and 3. 
However, if randomizing the time series of Wi{t)/W{t), which corresponds to 
eliminate all the possible correlations but keep the distribution form of step 
sizes, and then calculating the distribution of R{t) due to Eq. (7), one obtains 
the completely different results, that is, the convergence to Gaussian is faster 
and only the power law regions with Levy-like exponent a < 2 can be found, 
similar to the process of ordinary truncated Levy flight [26] . 

To directly calculate the cross-correlations between the relative updated wealths 
Wi{t)/W{t), we use the autocorrelation function of 

rr^^(T\ - {^{t)x{t + T))-{x{t)){x{t + T)) 



for some variable x{t). Although there is no correlation of relative wealths 
{x = Wi/W in (8)), the nonzero and positive correlation of the square {wi/WY 
is found with long persistence, as shown in Fig. 4, verifying the existence of the 
coupling between step sizes of the walk (7) for this multiplicative system. These 
correlations are attributed to the social term avb introduced in our system (6), 
supplying the coupling between different wealths. A larger value of a leads to 
an enhancement of cross-correlations, and then of a > 2 behavior, including 
the value of a and the extension of power- law region (see Fig. 3). Even when 
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Fig. 4. Correlations of {wi/W)'^ and R'^{t = 1) for N = 500, a = 0.0002, and A in 
the range (0.95,1.05) (averaged over 1100 runs). Inset: Correlations of (wi/W)'^ (av- 
eraged over 1000 runs) for N = 500 and different a values and A ranges: a = 0.0008 
and A G (0.9, 1.1) (circles), a = 0.0002 and A E (0.95, 1.05) (stars), and a = 5 x 10"^ 
and A G (0.98, 1.02) (triangles), which all correspond to the same exponent = 1-5 
of wealth distribution, as well as very small a = 2 x 10^^ and A € (0.995, 1.005) 
(diamonds). 

keeping the value of for wealth, the correlation of {wi/W)"^ increases with 
larger a. In the inset of Fig. 4, different values of coupling parameter a and 
different range of A are chosen so that the corresponding a^j keeps the same 
value 1.5, and from bottom to top the correlations increase with the increasing 
values of a. 

As the consequence of cross-correlations between Wi/W, the volatility cor- 
relations for R{t) (which is the weight-sum of Wi/W shown in Eq. (7) and 
equivalent to return r(r)), defined as the autocorrelation of square of returns, 
i.e., X = R^{t) in Eq. (8), exhibit the long-range persistence, as shown in Fig. 
4 for r = 1 (lower curve). Fig. 5 gives the results of volatility correlations for 
T = N = 500, which exhibits the slow decaying and seems to be exponential- 
type (inset of Fig. 5) rather than the power law decay found in real market 
data [28]. For correlation of return R the value is around 0, consistent with 
empirical observations [12,15,17]. This property of volatility clustering can be 
also reproduced in some microscopic models, but there an extra dynamics like 
diffusion of traders or feedback between price change and trading activity is 
needed [29], while in this asynchronous system it is obtained intrinsically. 
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Fig. 5. Autocorrelations of i?^(r = N) (volatility correlations, averaged over 1500 
runs), for = 500, a = 0.0002, and the uniform distribution of A in the range 
(0.95,1.05). The autocorrelation of R{t), not shown, are on the scale of this figure 
just a fiat line at high zero. Inset: Semi-log plot of volatility correlation. 

5 Discussions and conclusions 



It is interesting to relate the scales r and N here with that of the real markets. 
The more natural time unit of measurement is r = iV, which corresponds to 
average updating interval of stocks or agents. To describe the properties of 
market index, the relevant degrees of freedom could be interpreted as the 
major companies capitalization in the market, while for individual stocks A^ 
may correspond to the important agents dominating the stock price. Thus, 
the somewhat different behavior (value of a, exponential cut-off, etc.) found 
in different markets and stocks may be attributed to different effective A^. 
Moreover, according to our results in Sec. 3 the effective A^ is not very large 
(that is, the finite size effect is significant) in real market, indicating a phe- 
nomenon that is also found in some microscopic models [30]: It is a limited 
number of important players, not millions of small traders, who determine the 
market. 

In summary, we have shown that the stochastic multiplicative model (6) repro- 
duces two different power law behaviors found in reality for the wealth distri- 
bution (2) and the return distribution (4), as well as the volatility clustering. 
The dynamics of returns is described by a walk with steps of sizes obeying a 
truncated Levy-like distribution, and in particular, having cross correlations. 
These cross-correlations between relative updated wealths Wi{t)/W{t) are ex- 
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pected to be the origin of the a > 2 behavior in the tail distribution and 
the long-range volatility correlations of returns, and can be attributed to the 
coupling term in multiplicative system (6). 
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